Field-induced superconductivity with enhanced and tunable paramagnetic limit 
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We demonstrate that in a superconducting multilayered system with alternating interlayer cou- 
pling a new type of nonuniform superconducting state can be realized under in-plane magnetic 
field. The Zeeman effect in this state is compensated by the energy splitting between bonding and 
antibonding levels. Such compensation mechanism at low temperature leads to the field-induced 
superconductivity. We discuss the conditions for the experimental observation of the predicted 
phenomena. 



There are two mechanisms of the superconductivity 
destruction by a magnetic field: orbital and paramag- 
netic effects [J, Usually it is the orbital effect that 
is more restrictive. However in the systems with large 
effective mass of electrons 0, or in low-dimensional 
compounds, like quasi-one-dimensional or layered super- 
conductors under in- plane magnetic field p|, the orbital 
magnetism is weakened and it is the paramagnetic effect 
which is responsible for the superconductivity destruc- 
tion. The Chandrasekhar-Clogston paramagnetic limit 
0, is achieved when the energy of the polarization of 
the normal electron gas, — x„£f 2 /2, equals the supercon- 
ducting condensation energy — iV(0)Ag/2 , where N(0) 
is the density of states of the normal electron gas, x n i ts 
spin susceptibility and Ao = X.76T C is the zero temper- 
ature superconducting gap. This gives the critical field 
H p = Ao/(v2/-*b) of the first order transition at T = 0, 
fi B being the Bohr magneton. Later Larkin and Ovchin- 
nikov 8] and Fulde and Ferrell (FFLO) predicted 
the appearance at low temperature of a nonuniform 
superconducting state with the zero temperature criti- 
cal field H^ B LO = 0.755A /^ S , i.e. somewhat higher 
than the paramagnetic limit H p . This prediction was 
made for three-dimensional superconductors. In quasi- 
two-dimensional superconductors the critical field of the 
FFLO state is even higher, namely H^ LO — Aq/h b 
|lO| . while in quasi-one-dimensional systems there is no 
paramagnetic limit at all The appearance of the 

modulated FFLO state is related to the pairing of elec- 
trons with opposite spins which do not have the opposite 
momenta anymore due to the Zeeman splitting. 

In this Letter, we demonstrate that in a ballistic super- 
conducting bilayer at low temperature and strong enough 
coupling t\ 3> Ao between the conducting planes, the 
paramagnetic limit is enhanced up to H c ~ t\/n B far 
above the usual limit H^ B LO = A /(j, B . More pre- 
cisely, a very unusual superconducting phase is settled 
between a lower and an upper critical fields given by 
H B H C = t\ ± HqHq jt\ and below a maximal temperature 
of the order of T^/ti. Thus one obtains field induced 
superconductivity above the lower field while the upper 
one provides the enhanced paramagnetic limit which may 



be tuned by varying the electronic coupling t%. This is 
due to the compensation of the the Zeeman splitting by 
the energy splitting t% between bonding and antibonding 
electronic states of the bilayer. As another important 
feature of this new phase, adjacent layers support op- 
posite signs of the order parameter. Note that such a 
so-called ir phase was predicted before |12l Il3| for the 
superconductor-ferromagnet multilayered systems where 
the atomic superconducting and ferromagnetic layers al- 
ternated and were weakly coupled in contrast to the 
present system. 
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FIG. 1: Multilayered system. The n' unit cell contains two 
superconducting planes <p n and ip n . The transfer integrals are 
very different, t% <C t%. The superconducting phase difference 
X between two adjacent planes can be either or tv. 

We consider a model multilayered system with a 
cristallographic structure similar to those of the high- 
T c superconductors 0], as shown in Fig. 1. Namely 
we assume that the electrons are confined in the 
atomic planes with the same zero-field dispersion relation 
£(p) = P 2 /2m — Ef, Ep being the Fermi energy. The 
transfert integrals between the planes are different t\ 3> 
t 2 , both being smaller than the Debye energy Hujd *C Ep 
. In this structure the n th unit cell contains two con- 
ducting planes, labelled by ip n and (p n . The coordinate 
in the plane is r and the cells are separated by a distance 
a along the z-axis which is chosen perpendicular to the 
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planes. We suppose that the Cooper pairing occurs in 
the planes and perform our analysis in the framework of 
the standard mean field BCS Hamiltonian, 



H = H^ +H v0 + H t + ^ J d 2 rA£ (r), 



(1) 



with 



H 4>0 = X^( p )^.« (P)^n,« (P) 



p,n 

l 



2 A nia0 (q)tp n<a (p + q)<0„,/3 (~p) + h.c, 
(P) (*lVn,c (P) + ^ (p))+fcc. 



where A is the BCS coupling constant, summation over 
repeated spin indexes a,f3 is implied, and A na p(q) is 
the Fourier transform of the superconducting order pa- 
rameter A njQ!| g(r) = A„(r)icr^a , a v being the second 
Pauli matrix. The operators ip n a (p) and ip n a (p) de- 
stroy one electron with spin a and momentum p, respec- 
tively in planes ip n and ip n . Our model includes both 
translationnal and jauge symmetry breaking. Indeed, the 
superconducting order parameter in the planes ip n and 
ip n are respectively given by A„(r) = Ae lq r+lKno and 
A„(r)e* x , where q and n are respectively the in-plane 
and the perpendicular modulation wave vectors. These 
vectors and the superconducting phase difference x must 
be determined from the minimum energy condition. 

The Gor'kov equations corresponding to the Hamil- 
tonian JJJ are solved exactly. In the small A limit, 
the linearized anomalous Gor'kov Green function F£ = 

tp\ (p)V4 ( — P)/ reduces to the form 



F\ 



|f/)K -\t k \ 2 ) 



f ~2 



(2) 



where t k = t\ + t2e lka , ui±(p) = iu> ± £(p) — [i B H, t — 
t k +K and u)±(p) = uj±(p + q), uo being the Matsubara 
frequencies. 

We first solve the isolated bilayer problem = 
0. After integration of Eq.JJJ over £(p), one ob- 
tains the anomalous Eilenberger propagator /^( v f) = 
J(de/ir)Ft=£ 0=± i/J,„with 

fla = u + ifi B H + t^a + iat ie -'x/2 cos( X /2) 
A 2(uj + ifi B H + i^p)(u> + iii B H + i^ +iah) 

(3) 

for positive cj, where a = ±1 labels bonding and an- 
tibonding states while Vf is the Fermi velocity vector 
in the plane. Then the self-consistency relation (@J be- 
low implies that the superconducting phase difference x 
between neighboring layers is either zero or tt. In the 



absence of Zeeman splitting H = 0, the superconduct- 
ing order parameter is naturally the same in both layers, 
namely x = 0. In this case, /t = A /id coincides with 
the well-known anomalous Green function of a supercon- 
ductor in the limit A — > 0, and the self-consistency rela- 
tion gives the bare critical temperature T c of the isolated 
layer. More generally for x = and finite H, the in- 
terlayer coupling t\ drops from Eq.Q and one retrieves 
/<£ = A/(uo + ifj, B H) for a two dimensional supercon- 
ductor under parallel magnetic field. The other possible 
choice is x = tt, when the superconducting order param- 
eter is opposite on adjacent layers. For small values of 
the magnetic field, this later n phase exhibits naturally 
a lower critical temperature than the x — phase. 
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FIG. 2: Excitation spectrum. Usual singlet pairing (thin line 
circles) between opposite-spin electrons occupying the same 
orbital is affected by Zeeman effect. In contrast, n coupling 
(thick line) between two electrons occupying a bonding and 
an antibonding orbitals may lead to the cancellation of the 
Zeeman splitting. 

However for relatively large interlayer coupling t\ > T c 
and high field, the situation becomes drastically different. 
Indeed, the excitation spectrum consists of four different 
branches e = ±£ + [i B H ± t-y in the limit A — > 0, see 
Fig. 2. The singlet pairing may occur here between one 
electron in the bonding orbital and the other electron in 
the antibonding orbital. This results in a very special 
coupling where, if fi B H — t%, the Zeeman splitting is 
exactly compensated. Therefore enhanced superconduc- 
tivity is expected in the vicinity of n B H = t\, at least at 
zero temperature. 

In order to derive rigourously this prediction, we anal- 
yse the self-consistency relation, 



A = 2^|A|A^(0)t]T Re(/t >a (v F )), 



(4) 



w>0,a 



in the tt phase where /^(vj?) 
t\ in the following way 



/L(vf) = 



depends on the coupling 



(5) 



2(lu + i[i B H + iati + ivp ■ q/2) 



Here N2d(0) = to/ (2ir) is the two-dimensional density of 
states per unit surface and per one spin orientation, and 
the brackets (...) denotes averaging over the polar angle 
#=(vF,q). 
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We first discuss the zero temperature second-order 
phase transition between the normal metal and the 
7r phase, as a function of the magnetic field. From 
Eas. l|4l5|) . the critical field H is shown to satisfy 



^ B H/(2^T C ) 



Am 



(6) 



where X = |q| Ui?/(2^ s ), and Hq = Ao/2^ B is the 
critical field of the second-order superconducting phase 
transition in a two-dimensional monolayer. One must 
then find the value of X which maximizes the critical 
field H. If the 7r phase is assumed to be uniform in- 
side each plane, namely if q = 0, Eq.© merely re- 
duces to \H 2 — t\/ii 2 B \ = Hq and we obtain a lower and 
an upper critical fields respectively given by [i B H c = 
ti ± [i B HQ/2ti, in the limit t\ fJ, B H - Thus at 
zero temperature and strong enough coupling t\ 3> A , 
the superconductivity destruction follows a very special 
scenario. At low fields, superconductivity is first sup- 
pressed in the usual manner at the paramagnetic limit 
^2D L ° ~ Ao/m_b leading to the normal metal phase. 
Then further increase of the field leads to a normal to 
superconducting phase transition at the lower critical 
field. This superconducting tt phase is finally suppressed 
at the upper critical field. This is a new paramagnetic 
limit which may be tuned far above the usual one merely 
by choosing the coupling t\ greater than Ao. Thorough 
analysis of Eq.© shows that the upper critical field is 
even increased by an in-plane modulation in analogy with 
the two-dimensional FFLO phase [Io|. In this sense, the 
low-temperature tt phase in the present case should be 
called a FFLO-7T phase. The upper critical field is max- 



imal for the choice X = 
then Eq.® reduces to 



|q| v F /2fi B = \H-ti/fi B \, and 



\H - h/» B \ . H + h/vtg + 2y/Hh/^ B = 4iJ 2 , (7) 

which gives the upper and lower fields /J. B H C = tx ± 
(igHg/ti in the t\ 3> A limit. Note that the period 
of the modulated order parameter |q| _1 = £ (£i/A ) 
is larger than the corresponding period in the two- 
dimensional FFLO phase which coincides with the bal- 
listic coherence length £ = vp /Ao [Io| . 

Furthermore one may derive the full temperature-field 
phase diagram using Eas. (|4l5|) and the result is shown in 
Fig. 3. When the temperature is increased, the lower 
critical field increases whereas the upper one decreases. 
Along the upper (resp. lower) critical line the FFLO 
modulation is lost at some temperature T* p (resp. T* ow ). 
For higher temperatures a uniform n phase (U-7r) is re- 
covered and the temperature dependence of the critical 
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FIG. 3: Phase diagram for ti = 3A . U-x (resp. FFLO- 
X) denotes the uniform (resp. modulated) superconducting 
state with phase difference x between the planes. Thick (resp. 
thin) solid lines represents second-order transition between U- 
X (resp. FFLO-x) and normal metal phase (N). We expect 
the U-x/ FFLO-x transition lines (not calculated) to be in 
the vicinity of the (virtual) first order U-x/normal metal lines 
(dash-dotted). 



field is given by 
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1 . fi B H c (T) + ah 

2 + * 2ttT 



(8) 

where "J (x) is the Digamma function and ^(1/2) = 
-C - 2 In 2 ~ -1.963, C being the Euler constant. Fi- 
nally the lower and upper critical lines merge at field 
H c = t\jn B and temperature Tm — Ke~ c T 2 / (Ati) in 
the limit t\ 3> T. Therefore the field induced 7r super- 
conductivity is confined to temperatures lower than Tm- 
The structure of these U-7T and the FFLO-7T phases is 
reminiscent of the corresponding U-0 and the FFLO-0 
phases although the former are shifted to higher fields 
and lower temperatures than the later. 
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FIG. 4: Phase diagram for t\ = Ao. Solid and dash-dotted 
lines have the same meaning than in Fig. 3. 

Above results were obtained for relatively strong cou- 
pling. For lower coupling t\ ~ Ao, the U-7r and the 
FFLO-7T phases merge continuously into the usual \ = 
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phases as shown in Fig. 4, and finally disappear for t\ 
slightly smaller than Ao- From an experimental point of 
view, one might choose a system with intermediate cou- 
pling t\ small enough to settle the 7r phase island in an 
available range of temperatures but also large enough to 
separate the ir phase island from the usual superconduct- 
ing phases with x = 0. 

Hereafter we discuss various physical mechanisms lim- 
iting the above predicted ir superconductivity: role of 
finite t<i, impurity and orbital effects. 

For finite coupling t% between the bilayers, the self- 
consistency relation (@J together with Eq. (0 leads to the 
following equation 



hi 



£ ^(H B H + X cos 0)2 - (t[ + at 2 f 

a=±l 



0. (9) 



For X = \H — (ti +t 2 )/fJ. B \, the singularity around 
9 = t: produces a corrective term to the upper and lower 
critical fields found previously. However this correction is 
negligible if tit 2 <C Mb-^o- ^ n tne opposite regime, that 
is for larger values of the inter-bilayers coupling t 2 , the 
bonding and the antibonding electronic levels form bands 
whose dispersion avoids exact compensation between the 
intra-bilayer coupling t\ and the Zeeman splitting. When 
<i = <2, the quasi- two-dimensional case is retrieved: 
the presently studied tt phases are lost in favor of FFLO 
phases modulated either along the planes or perpendic- 
ular to the planes. 

Hence we see that the interlayer coupling must be 
rather small to prevent the smearing out of the elec- 
tronic levels participating to the tx coupling. Extremely 
low coupling may in principle leads to the suppression 
of the transition due to the two-dimensional fluctuations 
[lH. However in the limit t 2 > T c yjT c /E F the fluctu- 
ations are limited in a very narrow temperature region 
near the critical temperature 01- For t 2 <C T C ^T C JE F 
Kosterlitz-Thouless regime may be attained near T c |l8| 
but the long range order is restored outside the vicinity 
of the critical temperature, namely for (T c — T)/T c > 
(T C /E F ) ln(T c 3 / 'E F t 2 ). This weak logarithmic divergence 
of the fluctuations as t 2 — > means that in practice very 
weak coupling is enough to restore the transition. Finally, 
the FFLO-7T phase is fully established if t 2 <C A 2 ,/^. 

Impurities produce further broadening of the electronic 
levels over an energy range 1/r, r being the elastic colli- 
sion time. One may infer that disorder leads to the de- 
struction of field-induced superconductivity since the 7r 
phases originate from exact compensation of the Zeeman 
splitting by the t\ splitting between bonding and anti- 
bonding electronic levels. As a rough estimate, 1/r plays 
a similar role than the inter-bilayer coupling t 2 and thus 
1/r = fj/gHg/ti defines a degree of disorder above which 
the 7r phases are unlikely to survive. Consequently the 



observation of the above predicted phenomena requires 
even cleaner samples than the observation of the two- 
dimensional FFLO state which was expected to exist for 
1/r < n B H 0. 

Up to now the orbital effect has been neglected. If we 
choose the in-plane coherence length as £ « 100 A and 
ti/Ep ~ 10~ 2 — 10~ 3 , orbital effects are expected to be 
important above H° c ^ b = (</> /£ 2 )(£ F /h) « 10 - 100 T, 
4>q being the magnetic flux quantum. 

Modern nanotechnology permits now to design and 
fabricate a rich variety of different superconducting struc- 
tures. In the present Letter we demonstrate that the sys- 
tem comprising weakly coupled superconducting bilayers 
may reveal a new type of the FFLO-7T superconducting 
state when the spin effect prevails against the orbital ef- 
fect of the magnetic field. The critical field of such state 
can exceed many times the standard paramagnetic limit. 
Moreover the phenomenon of the field-induced supercon- 
ductivity at low temperature is predicted. Note that 
this field-induced superconductivity is very different from 
that recently observed in quasi-two-dimensional organic 
conductor A — (BETS) 2 FeCl4 0, where the exchange 
field of aligned Fe 3+ spins compensates the external field 
by means of the Jaccarino-Peter effect • 

We are grateful to M. Houzet and A. Koshelev for use- 
ful comments. 
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